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Abstract. We introduce a designer approach for extended Brans-Dicke gravity that allows 
us to obtain the evolution of the scalar field by fixing the Hubble parameter to that of a 
uiCDM model. We obtain analytical approximations for 0 as a function of the scale factor 
and use these to build expressions for the effective Newton’s constant at the background and 
at the linear level and the slip between the perturbed Newtonian potentials. By doing so, we 
are able to explore their dependence on the fundamental parameters of the theory. 
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1 Introduction 

Over the next decade, we expect a step change in our understanding of gravity on cosmological 
scales. Surveys of large scale structure should be able to pin down the expansion of the 
Universe and the growth of structure with exquisite precision [1-4]. These new data sets 
should allow us to constrain modifications to general relativity at a level which may be 
comparable to those obtained on astrophysical scales. 

If we are to fully take advantage of these data sets, it is essential to have a detailed 
and accurate understanding of how different observables depend on our assumptions about 
gravity. In particular, we should know how deviations from general relativity will affect our 
observations: whether the effects are large or small (given what we know on astrophysical 
scales) and how correlations between the observables themselves might be indicative of some 
underlying structure. 

There has been a formidable campaign to develop methods for studying the effects of 
modified gravity on large scales (for a compendium of theories, see [5]). A different approach 
has been to develop a unified method of parameterizing all possible theories at the linearized 
level (for a selection of methods, see [6-9]). Yet, while there is an inexorable momentum 
that should lead to a battery of effective techniques for extracting useful information from 
the data, we do not have yet a firm understanding of what to expect. By this we mean that, 
given certain theoretical assumptions, what our observables should look like, i.e. what values 
should they take and how should they be interrelated as a function of whatever fundamental 
parameters we might consider. 

In principle, the step from taking the parameters, (with i = 1, • • • N ), of some underly¬ 
ing theory and working out the resulting phenomenological parameters, /3j (with j = 1, • • • M ) 
tied to observations, should be straightforward. In practice, the process can be complicated, 
highly non-linear, degenerate and normally obscures the relationship between the prior as¬ 
sumptions on ai and the resulting theoretical priors on /3j. One way around this is to develop 
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an approximate mapping between the two sets of parameters and, wherever possible, analytic 
relations between the two. Furthermore, if one can find a method for restricting the range of 
a* given some assumptions about a subset of the /%, one can quickly surmise what correla¬ 
tions and covariance one should expect for the remaining phenomenological parameters. In 
this paper we propose an approach to do so, considering a restricted model for cosmological 
modifications to gravity. 

Our starting point is a well known theory, the Brans-Dicke (BD) theory of gravitation 
[10]. This theory is the simplest scalar-tensor theory one can envisage [11-15] and is consid¬ 
ered a viable alternative to General Relativity, one which respects Mach’s Principle. Since 
its formulation, this theory has been exhaustively studied as a possible alternative solution 
for the accelerated expansion of the Universe. 

It has been shown that Brans-Dicke theory can produce accelerated solutions for small, 
negative values of the BD parameter cubd [16, 17]. Given that one recovers standard GR 
in the limit where cubd —> oo, such values of the cubd clash with Solar system constraints 
[18, 19]; furthermore, recent constraints with the latest CMB data are also not compatible 
with such low values of cubd [20, 21]. Several modifications of this theory try to include 
self-interacting potentials [22-24] or consider a field-dependent Brans-Dicke parameter u(4>) 
[25], without solving this problem successfully. Also, models with a non-minimal coupling of 
the scalar held have been considered in Refs. [26-29]. 

In this paper we construct a theory of designer , extended Brans-Dicke gravity and use 
it to characterize the form of the observables we might measure. This theory is ”extended” 
because we include a potential for the Brans-Dicke held and we dub it ’’designer” (the term 
“designer” was hrst used in models of inflation that attempted to match observations by 
designing the density fluctuation spectra [30]) because we reconstruct the potential (which 
might not have an analytic form) from a desired background evolution. While such a theory 
does not seem fundamental, it might be seen as an approximation to a scalar-tensor the¬ 
ory which has a particular, a priori, form of the background evolution. Our construction 
allows us to hnd a number of analytic approximations and, in doing so, lets us gain a firmer 
understanding of the phenomena we want to study. 

Our designer approach for the extended Brans-Dicke gravity is novel. It allows us to 
retrieve the evolution of the scalar held, (j>, by hxing the background evolution and is robust 
for high values of the BD parameter, which is the regime we are interested in. This method 
not only works for a ACDM like evolution with an effective equation of state = — 1, but 
is also applicable for models with w e Q > —1 as in a teCDM scenario. And, for both cases, we 
are able to retrieve analytical approximations for <j> as a function of the scale factor a which 
could prove useful for a faster and more efficient htting of models to data. 

The paper is structured as follows. In Sec. 2 we introduce the Brans-Dicke theory with 
a constant wbd parameter. In Sec. 3 we describe the designer approach, motivated by an 
analysis of the behavior of this theory when we have a constant potential V((pi). In Sec. 4 find 
approximate analytic solutions to the evolution of the scalar field and use it to infer the shape 
of the potential. We then use these results in Sec. 5 to construct analytical approximations 
to the phenomenological parameter which can be constrained by data. In Sec. 6 we discuss 
our results. 
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2 Extended Brans-Dicke gravity: background equations 

The action for extended Brans-Dicke theory in the Jordan frame, is given by 

S = ^2 j d^x\/—g (d>R - ^ (<90) 2 - 2V(4>)^j + S m , (2.1) 

where S m [$ m ; g^] is the minimally coupled matter Lagrangian and k 2 = 8irG, where G 
is Newton’s gravitational constant measured today. Varying the action with respect to the 
metric elements, we find the Einstein equations, 


K rpm , W BD 


r< — rpm , 

^ (IV — , ± my I 




4 > ,fJ.4 > ,v T , dfiuGicj)] 9pv, 


( 2 . 2 ) 


where T™ u is the matter stress-energy tensor. 

By varying the action (2.1) with respect to the held, one gets the field’s equation of 
motion 

a * = rpk~ B -iT < 2 - 3 > 

where V0 = dV/dcj). Considering a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) met¬ 
ric, ds 2 = — dt 2 + a 2 (t)dx 2 , this equation reads 


■ + 3 Hcj) = 


K 2 p m 4V(0)-20^ 


3 + 2wbd 


+ 


3 + 2wbd 


(2.4) 


where p m is the matter’s energy density and H = a/a is the Hubble parameter. The latter 
is determined by the two Friedmann equations, which are written as 


3H 2 (f> = K 2 p m -3 H<j ) +‘^U + V M 

2 (p 

2 H + 3 H 2 = -k 2P ^ _ _ 2 + 

(b 2 6 Z <p (b <b 


(2.5) 


Lastly, from the previous equations, one can define an effective equation of state for the 
dark energy component of our model, which is given by 


w e s = 


0 2 w(0) + AHij) + 20 - 2V(0) 
0 2 w(0) - 6H<j) + 2V(0) : 


( 2 . 6 ) 


where w(0) = cubd/ 0 and, even more straightforwardly, one can define the fractional effective 
dark energy density parameter, 

P<t> 


where the effective energy density is given by 

wbd 0 2 


3 H 2 cb' 


P<t> — 


- 3FT0 + V (</>). 


(2.7) 


( 2 . 8 ) 


- 3 - 

















2.1 Constant Potential V(0) 

Before proceeding to the designer approach, we can get an idea of the different effects at 
play in extended Brans-Dicke gravity by considering the case of a constant potential U(0). 
For all our calculations in this section, we have U(0) = 3H q (1 — fi m ) = V, where D m is 
the fractional present-day energy density of matter. For a perfect ACDM scenario we should 
have an effective dark energy equation of state equal to —1 during the whole cosmological 
evolution, with the scalar held remaining perfectly still and showing no evolution at all. 
However, in the Brans-Dicke paradigm, the held should always evolve even if its dynamics 
are subdominant (in “slow roll”) compared to the potential V. Hence, effectively, we will 
have a quasi-ACDM evolution. 

We start by numerically solving the scalar held evolution using Eqs. (2.4) and (2.5) 
considering a constant potential as dehned in the previous paragraph. We set the initial 
conditions for the scalar held deep within the matter dominated regime at a redshift around 
Zj ~ 1000. For this, we consider a known solution of Brans-Dicke gravity given by [31-33] 

0 = 0oa 1/{u} BD+1) , (2.9) 


where 0o = (2wbd + 4) / (2ljbd + 3). This solution is, in fact, an attractor solution of the 
system derived in the absence of a potential U(0) and for a Universe dominated by matter 
alone [31-33]. The scale factor, on the other hand, evolves as [31-33] 


a{t) 


^ \ (2oj B d+2)/(3oj B d+4) 

to) 


( 2 . 10 ) 


and we see that, in the GR limit of wbd — > oo, 0=1 and a(t ) oc t 2 / 3 throughout the 
matter dominated regime; to is related to the inverse of the present-day value of the Hubble 
parameter, Hq, such that fo-^o = (2 wbd + 2) / (3ljbd + 4). The value of 0o ensures that, in 
a matter dominated Universe, we would measure an effective gravitational constant today, 
G e g, equal to the actual Newton’s gravitational constant, G, in Cavendish-like experiments. 
This assumes, of course, that the Solar system value of 0 is representative of the Universe as 
a whole, which may not be entirely accurate [34]. 

Let us also point out that, in a matter dominated flat Universe, the matter density 
will not be precisely equal to the critical density due to a very small, negative, and almost 
negligible contribution from the scalar held dynamics. It is possible to rescale the matter 
density (as in Ref. [35]), but we opt not to do so, since the correction is negligible in the 
cjbd >> 1 regime we are mostly interested in this work. 

In Fig. 1, we have the numerical evolution of the scalar held plotted against the power- 
law solution given by Eq. (2.9). We can clearly observe that, even in the presence of a constant 
potential V, the Brans-Dicke scalar held evolves according Eq. (2.9) at early-times, during 
the matter dominated epoch. Only at late-times, close to a = 1, we see a slight departure 
from the power-law of Eq. (2.9), when the dark energy component begins to dominate and 
accelerates the scalar held. 

Still in Fig. 1 we can observe the numerical evolution of the dark energy effective equa¬ 
tion of state iCeff as given by Eq. (2.6). We observe a very sharp transition from —0.4 to —1 
that we will explain later on. For now, we can conclude that, even though the scalar held is 
accelerated by the presence of the constant potential V, its dynamics remain subdominant 
(the aforementioned slow roll evolution) and allow for a late-time potential dominated epoch 
with = — 1 
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Figure 1: We plot the numerical evolution (solid lines) of the scalar field evolution and the 
effective dark energy equation of state, w e s in the presence of a constant potential. We note 
that, in the matter dominated regime, (j) evolves according to a known power-law solution 
given by Eq. (2.9), which we also plot (dashed lines). 


Having shown in Fig. 1 that we recover the power-law solution given Eq. (2.9) at early- 
tinres, we now extend its application by using it in the effective equation of state rc e ff given 
by Eq. (2.6) in the presence of a constant potential V. Hence, we approximately obtain 

_ 4 - 4u;bd Va 3 /H$ 

WeS ~ —10 + 4u; BD Ea 3 /’ 

in the limit of wbd >> 1, and where we have also used Eq. (2.10). Hence, in the matter 
dominated regime, the potential contribution is suppressed by the scale factor leading to 
w e fj ~ —0.4 (unless cubd -> oo and V / 0). Thus, for values of wbd which are consistent 
with Solar System constraints, it is impossible to get an accelerated solution without adding 
a potential V((/)), that may not necessarily be constant. However, with a constant potential 
V(cf>), one gets w e g = —1 at late times after a sharp, non-smooth transition from w e s ~ —0.4, 
which we have seen in Fig. 1. 

An effective equation of state ~ —0.4 at early times could constitute a problem, 
eventually compromising the extension of the matter dominated regime and rendering the 
model inviable. However, calculating given by Eq. (2.7), explicitly during the matter 
dominated regime using Eqs. (2.9) and (2.10), one gets 
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which, for large values of wbd is negligible at early times. 

Also, we note that the discontinuity in w e s happens due to a zero crossing of the 
denominator of Eq. (2.6). If we change from physical time t to the natural logarithm of 
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the scale factor, dt — > din a, we have that d/dt — > Hd/dlna. Therefore, neglecting the 
cj)' 2 (the prime denotes a derivative with respect to In a) term because this is proportional 
to (1 + <ubd) - 2 before the transition and for large wbd, the denominator of w e g- can be 
approximated to just — 3c/) 1 + V(4>)/H 2 . Therefore, given that, in the matter dominated 
regime, V{(j))/H 2 oc V(cj))a 3 / Hq is an increasing function of the scale factor, there will come 
a point at which this term will be equal to 3<p', leading to the discontinuity in rc e fj. For the 
constant potential, the scale factor of the discontinuity is approximately 


f »m 1 y /3 

V 1 — fl'm 1 + WBD J 


(2.13) 


The discontinuity in w e g has no impact on the background expansion of the model: if we 
take the second Friedmann equation and p m = 0, we have 

- = (l + 3 WeS 

O 2 \ Pm/ P<j> T1 

Since the divergence in w e g happens due to p<f, crossing zero, as we just discussed, no diver¬ 
gence is seen in the evolution of d because the term p m /p<p follows the behavior of ru e ff. 

Finally, only when wbd — > oo (the General Relativity limit) does one get tt> e fr = — 1 
throughout the whole evolution, as seen in Fig. 1. Here the potential E(</>) will dominate and 
the scalar field dynamics is heavily suppressed. The discontinuity in ru e ff will now happen at 
a much earlier time, as is clear from Eq. (2.13), leading to a smooth u> e g = — 1 in the case of 
a constant potential. 


(2.14) 


3 Designer extended Brans-Dicke gravity 

Having presented the general form for extended Brans Dicke gravity, we now proceed to 
construct an algorithm that will lead to a particular expansion rate or, more specifically, 
to an effective equation of state. Hence, effectively, we design and impose the background 
history we wish for our model which in turn determines the dynamical evolution of the Brans- 
Dicke scalar field. We note that the authors of Ref. [36] suggested the designer method we 
will describe further, but did not fully explore its consequences. 

Following the previous section, we have shown that, at early-times, the scalar field 
will follow the matter domination attractor solution irrespective of the presence of a scalar 
potential V(4>). At late-times, its evolution should be dominated by leading to a 

departure from the matter dominated attractor solution. Therefore, we now try fixing the 
background evolution to match that of a standard flat tcCDM scenario, such that 

H\a) = ^ [H m a- 3 + E eS (a)] , (3.1) 

where Q m is the present-day fractional matter energy density, and the dark energy component 
will be fixed as 

E eS (a) = (1 - n m ) e 3/ a 1 (i+^eff)dina_ (3.2) 

We will be assuming that the effective dark energy equation of state u; e ff is a constant such 
that w e g > —1. We should be clear, however, that this is not a limitation of this procedure: 
it can be easily extended to a varying u; e ff by providing a as a function of the scale 
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Figure 2: We show the evolution of <j) and w; e fr as a function of the scale factor, for 
wbd = 1000. On the left, the top plot has w e s = —0.80 today, while the bottom plot has 
= 0.30 and different values of w e g today. On the right we show the respective evolution 
of rCeff. 


factor a. We merely choose to do so in hope of finding analytic expressions for some of 
the observables in terms of the fundamental parameters of the theory. Therefore, we can 
now numerically evolve the scalar field just by using Eq. (2.4) without evolving the Hubble 
parameter using Eq. (2.5). We are also effectively parameterizing Eq. (2.8) so that our dark 
energy component’s energy density matches a tcCDM type and are not worried with its exact 
numerical evolution. 

We then take the approximation of considering the scalar field potential to be determined 

by, 

V{<t>) = 3 Hi (1 - H m ) e 3 / Q 1 (i+^ef I )dina ) (3.3) 

meaning that we are considering that the main contribution to the effective dark energy den¬ 
sity comes from the scalar field potential, with the scalar field dynamics being sub-dominant. 
With this approximation we also don’t expect to affect the matter domination attractor so¬ 
lution at early times since, as seen before, the potential contribution to is not relevant in 
the matter dominated regime. 

To generate our numerical results we have fixed the initial value of the scalar field <j>(zi) 
and < t>'(zi ) to match the matter dominated attractor solution value at a redshift of Zi = 1000. 
In Fig. 2 we plot the evolution of 4> and w e s for different values of Q m , w e s, and cubd by 
numerically solving Eq. (2.4) and fixing the evolution of H with Eq. (3.1). 

We note that, similarly to what we observed in the constant potential case, the presence 
of the dark energy component leads to a departure of 4> from the matter domination attractor 
solution at late times, leading to a scalar field value higher than cf >o at the present. And Fig. 2 
makes it clear that this departure happens earlier in time and is more significant the earlier 
the dark energy component starts to dominate at late-times (which happens the bigger w e s 
is or the smaller is). This means that, the higher w e g is, the more relevant the scalar 
field dynamics becomes. Hence, our designer approach breaks down if w e g is much higher 
than —1. 
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Looking at Eq. (2.8), one might be concerned about the numerical evolution of the 
effective dark energy density which we parameterized by Eq. (3.2); we would probably not 
recover a flat cosmology today due to the contribution of the scalar field dynamics to the 
overall critical density of the Universe. If we were to compute p^ numerically with Eq. (2.8), 
one could adjust the weight of the potential V{4>) to compensate for the dynamics of the 
scalar field and recover = 1 — f2 m today. Hence, in effect, V((f>) = 3HQ$l ( f ) a-^ lJrWeS \ 
where 0,$ ^ (1 — H m ) could be found by performing a simple binary search, for example. We 
will provide an approximation for this factor using our analytical solutions for (f> in Appendix 
A. 


We can also study the evolution of rc e ff in Fig. 2. We see that we again have a sharp 
transition from the matter domination attractor regime w e s = —0.4 value at early times to 
the value we fix w e g to at late-times. The scale factor at which this transition happens can 
be estimated from 


®disc 


/ Q m 1 \ 3u 'eff 

V 1 — Hm 1 + WBD J 


(3.4) 


making it clear that, the larger w e s is and the earlier our dark energy component becomes 
relevant, the earlier this transition happens. Also, even though we don’t show that explicitly, 
we recover the GR plus u>CDM limit when we take wbd —> oo, and w e fj should then be 
equal to the value we fix it to be throughout the whole evolution, since the discontinuity now 
happens earlier or may even be completely avoided. 


4 Analytical solutions for (f> 

With our designer approach in hand, we can now proceed to find analytical approximations 
to the scalar field evolution which, in turn, can be used to construct approximations to our 
observables. We first consider the ACDM-like case and then generalize to an arbitrary (but 
constant) effective equation of state w e s. 

4.1 W e fi = -1 

We start by expressing the scalar field equation of motion, given by Eq. (2.4), in terms of 
In a. We then simplify it by simultaneously neglecting the <j)" and (j)' 2 terms, yielding 

$ (y 1 ^m«~ 3 \ _ 4(1 - O m ) +q~ 3 fl m 

\ 2 1 —H m + o _3 H m y d (1 — fi m ) + a _3 H m ’ 

where d = (2cubd + 3). The solution for the scalar field will be a fully analytical expression, 
given by 

4>{a) = (t>{o,i)g{ai)~ l g{a ), (4.2) 

where 4>(ai) is the scalar field value at a high redshift Zi set by the matter dominated attractor 
solution, or can be fixed to be ^>o at ai = 1. The function g(a) is given by 

g(a) = ad (2a 3 (1 — H m ) + H m ) 3d (4.3) 

We show the evolution of 4> predicted by this solution in Fig. 3. It exhibits a tendency 
to overestimate the deviation from the matter domination attractor solution at late-times. 
However, its errors are small, specially when considering the considerable simplification we 
have found to the full numerical analysis of our designer approach. 
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Figure 3: We show the analytical solution for the scalar field, 0 so i, predicted by Eq. (4.2) (in 
solid lines) for different values of wbd and in the upper plot. We compare the analytical 
solutions to the numerical evolution of 0 (in dashed lines) predicted by our designer method, 
and show the relative error in the bottom plot. The errors are shown in %. 


4.2 w eS ^ - 1 

We now extend our analytical approximation for cosmologies with w e fj ^ — 1. In these 
circumstances, we expect the dark energy component to become relevant earlier, and hence 
produce larger deviations from the matter dominated attractor prediction. We will focus 
mainly in the late-time evolution of 4>, when the dark energy component comes to dominate. 
For that effect, we re-express Eq. (2.4) in terms of In a, and assuming = l/<f>' dV/dlna, 
we approximate it as 


12 ^ (1 - n m ) + i8cj) (1 - n m ) (1 + w eS ) 
d(l-n m + n m a 3w es) 


30 ' 2 , 


(4.4) 


where d = (2wbd + 3) and we have also neglected terms proportional to 0", q i>' 3 and (1 + 
iCeff)(1 — fl m ), the last two arising with the derivative of H. We have not included the 
matter driving term that dominates at early-times. Assuming that the driving term from 
the potential slope is much more significant than the V (0) one- which effectively means <f>' is 
much smaller than unity for large wbd- we take the square root of this equation and perform 
a Taylor expansion of the left-hand side, obtaining: 


6^' (i - n m ) 


: + \/l8(l - n m ) (1 + rc e ff)(/>-0 / v / 3rf(l - ^rn + £l m a 3Wea ) ~ 0 (4.5) 


■\/ is (i — n m ) (i + w e ff) 

With these approximations, the solution for this equation is given by 

0(a) = 0(a i )/(a)/(a i )~ 1 , 
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Figure 4: We show the late-time analytical solution (solid lines) of the scalar field predicted 
by Eq. (4.6), (j> so \, when w e s / —1 for different values of wbd and fi m . We compare it to the 
numerical solution (dashed lines) predicted by our designer approach in the top plot. We 
show the relative errors in the bottom plot, in %. The red lines are for wbd = 100, while the 
green and blue lines are for wbd = 1000. 


where (j>(ai) is the value of the scalar field at a desired scale factor dj. This can either be 
set to the matter dominated attractor solution at a redshift Zi ~ 10 or to at a = 1 if 
one wants to fix the present-day value of the scalar field to recover G e fr/G = 1 today. The 
function /(a) is given approximately by 


/(«) 


T + x' 


V6Vd(l + w e ff ) 3 / 2 

™ eff (- 2 + 3d <’ 1 + u ' eff )) 


x — 1 


(4.7) 


where x = \J 1 + a 3 ^ eff and we have neglected similar terms whose exponents were 

proportional to dr 1 . In Fig. 4 we compare the late-time evolution of 4> predicted by Eq. (4.6) 
with the numerical evolution found in our designer approach. We do so by fixing ^(a*) to the 
matter domination attractor solution at = 10 for all the cases. We see that this solution 
works better for larger wbd- Nevertheless, even if the agreement with the numerical solution 
is not perfect, the errors are small, and the overall form of (f> is excellent for such a simple 
approximation. 

We are now also in a position to reconstruct the effective form of the self-interaction 
potential V((f>) across the entire cosmological evolution. For that, we invert the solutions for 
4> to get the scale factor as a function of the scalar field. We use the field’s matter dominated 
attractor solution at early-times and our analytical approximation at late-times. Hence, the 
potential will be given by 
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Figure 5: We show the reconstructed functional form of V(4>). The red lines have cjbd = 100, 
while the blue line has cubd = 1000 . 


3flg(l-n m ) (0/^>o) 3 ( 1+U,eff )( 1+WBD ) 5 during matter domination 

~3(l+^ off ) 

u ’eff 


V(<j>) 


3Hq(i — n m ) 


(-2+3d(l+m eff ))m eff 
2 2 / 3 (l-n m ) 1/3 (0/c) 3^Vd(l+^ eff ) 3/2 
/ (— 2+3d(l+ Weff )> eff \ 2/3 

1 ~{ 4 >/ C ) ^Vdll+^eff) 372 J 


,1/3 


at late-times, 


(4.8) 

where c = <?;(a,;)/(a,;) , as defined in Eq. (4.6). We plot the late-time form of the potential 

V((j)) and V^/V in Fig. 5 for different values of wbDi Well and fl m . We can observe that 
V(4>) exhibits a simple form, as in a standard run-away potential, with the slope decreasing 
at higher values of 0 or, equivalently, close to the present. We see as well that V^/V takes 
significantly high, absolute values. This justifies our assumption in considering just the effect 
of the slope of the potential in the evolution of 4> when w e s 7 ^ — 1. Indeed, this is the term 
that will have the most effect on the scalar field dynamics, leading to a significant departure 
from the attractor solution at late-times for w e s > — 1 . 

In Fig. 5 we can also observe that the slope of the potential becomes more significant 
for higher values of wbd- This seem to contradict what we have seen in Fig. 4, where the 
scalar held dynamics seem to be more relevant, the smaller wbd is. However, the source 
terms for the evolution of <f> are suppressed by a factor proportional to Hence, for a 

larger value of ubd, the only way to have significant held dynamics at late-times, and hence 
induce a signihcant departure from the matter dominated attractor solution that produces a 
tCefr 7 ^ —1, is to have a very large source term. Finally, we can also see how, for larger rc e ff, 
we recover a more tilted potential: the more relevant we set our dark energy component to 
be, the more signihcant we expect the scalar held dynamics to be at late-times. 
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4.3 A global solution 

In the previous two sections, we presented solutions for the evolution of the scalar field that 
worked well for w e .g = — 1 and > — 1 separately. We will now propose an approximate 
global solution: 

</> g iobai(a) = 4>(ai)f{a)g(a)f(a i )~ 1 g(a i )~ 1 , (4.9) 

which is just the product of the solutions we previously found for u; e ff = — 1 and for u; e fr > — 1. 

Note that when = —1, we have that f(a n ) = /(a) = 1. Hence, ^> g iobai(a) will be the 
exact solution for the scalar field equation of motion under the assumptions we discussed in 
Sec. 4.1. When w e s > —1, we note that the main contribution will come from the /(a) and 
/(oj) terms; we already have seen in Section 4.2 how the scalar field dynamics are more signifi¬ 
cant when u; e H- > —1. Not only that, but we note that dg{a) / dhia oc d _1 , whereas df(a)/d In a 

produces terms proportional to (^dj • Hence, assuming dg(a)/dlna « d/(a)/dlna when 
> —1, (^global(a) w hl be a solution of Eq. (4.5). We will use this full solution in the fol¬ 
lowing sections for the phenomenological parameters, and show that it is indeed a good 
approximation for the overall behavior of the scalar field. 


5 A model for the phenomenological parameters. 


It has been shown that, at the level of the background and linear cosmological perturbation 
theory, it is possible to completely characterize any modified theory of gravity in terms of 
a handful of time dependent functions [ 6 ]. We proceed to do so with our designer extended 
Brans-Dicke gravity. We have already discussed two of our time dependent functions: the 
time varying Newton’s constant (at the level of the background), Go = G/<f> and the effective 
equation of state, u; e ff- 

For linear perturbations, following the notation of Ref. [37], we consider a perturbed 
metric about the FLRW background in the Newtonian gauge, 

ds 2 = —(1 + 2 ^>)dt 2 + a 2 (t)( 1 + 2^)Sijdx l dx\ (5.1) 


where 4' and $ are the scalar perturbations that we will refer to as Newtonian potentials 
and are decomposed as a series of Fourier modes of scale k (h/ Mpc). 

If we are interested in the the impact of matter perturbations on galaxy and weak lensing 
surveys, we can focus on the modes that are well within the Hubble radius, i.e. such that 
the condition k 2 /a 2 H 2 is respected. In this quasi-static regime the evolution equation for 
the matter density perturbation 5 m can be approximated as [38] 

C + { 2+I p)C- ~ 0, (5.2) 


where fl m (a) = Pm/3H 2 , and G e ff/G will be dependent on the model. The primes represent 
derivatives with respect to In a. In the extended Brans-Dicke theory, G e g- is given by [37] 


G e ff 14 + 2 cjbd + 2<^ (A/a/Zc) 2 
G (j) 3 + 2cjbd + 2</> (Ma/k) 2 

where the M term is [37] 


M 2 


Vc/xj) + 


WBD 

(j ) 3 



+ 3 


(5.3) 


(5.4) 
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At late times, when the dark energy component starts to become relevant, the mass term 
can be simplified and expressed in terms of the potential alone using the scalar field equation 
of motion, such that M 2 ~ [37]. where and correspond to the second and 

first order derivatives of the potential with respect to the scalar field, respectively. 

One can also define the gravitational slip 77 corresponding to the ratio between the two 
Newtonian potentials [37] 

£ _ 1 + ^bd + ct>(Ma/k) 2 

* 2 + wed + <t> (Ma/k) 2 


which, again, should depend on the specifics of the scalar-tensor model. Lastly, the sub¬ 
horizon version of the Poisson equation can be written as [37] 

k 2 

^ ~ —47rG'effPm<5 m - (5.6) 

a z 

In standard GR, when we neglect matter shear, the anisotropy equation between the 
Newtonian potentials becomes a simple constraint equation, T = ( I>, and 77 should be 1 
throughout the cosmological evolution, as should G e g/G. Hence, in a modified gravity theory, 
a deviation in these parameters signals a departure from standard GR that can potentially 
be measured. From Eqs. (5.3) and (5.5) it is clear that the GR limit is recovered when 
wbd —> oo, as expected, or when the field becomes supermassive and M 2 —> oo. But we now 
wish to understand how these functions depend on time. To do so, it is convenient to study 


£qs 

VQS 


lim 

k—> oo 


G e g 
G 


lim 77 = 

k —^00 


14 + 2wbd 
</> 3 + 2wbd 

1 + WBD 

2 + wbd 


(5.7) 


and the inverse length scale 

k M = \ r~~- - Ma ■ (5-8) 

V 1 + cubd 

From Eq. (5.7), we see that t]qs is constant throughout the cosmological evolution, 
independently of the scalar field dynamics [36-38]. Its GR limit is trivially recovered when 
we take cubd —> 00 . 

On the other hand, the late-time evolution of the mass scale parameter, kM, can be 
written as: 


k M (a) « ~^=H 0 ^2 (1 - H m ) (1 + w e ff) a + a 1 Q m (2 + w e g) (5.9) 

which is valid for w e g > — 1 and large wbd- In the limit w e g = — 1, this equation predicts a 
non-zero value for kM, whereas it should be exactly zero throughout, as predicted by explicitly 
using our global solution for (j) in Eq. (5.8). This should be evident as V^> = R/+ = 0 when 
w e g = — 1. This approximation works fairly well for small redshifts and better for larger wbd , 
as can be seen in Fig. 6 . 

We can also observe that k m is a fairly negligible quantity, corresponding to scales 
which are of order or greater than the cosmological scale. To do so we compare kM to 
the comoving horizon, aH = iLox/Hm ® -1 + (1 — H m )a _ ( 1 +3«' e fi). it is not hard to see that 
kM/(aH) < 1. Hence, the scale at which k/kM becomes relevant is approximately the 
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Figure 6: We show the evolution of Izm- The solid lines are the numerical evolution obtained 
using the global solution for cj) given by Eq. (4.9), while the dashed lines show the evolution 
predicted by the approximation given by Eq. (5.9). The red lines have wbd = 1000, the 
blue lines are for wbd = 10000 and the green lines have u>bd = 1 x 10 5 . We have fixed 
Hg = 1/2997.9 h/ Mpc, with h = 0.68. The black lines are for w e s = — 1, and the results 
shown are independent of cubd- 


same at which the perturbations fc-modes become sub-horizon, which is at the basis of our 
assumptions. Therefore, taking k/ku >> 1 is an excellent approximation on quasi-static 
scales. 

To understand the parameter dependence of £qs we perform a Taylor expansion around 
a = 1 using our approximate global solution for the scalar field. We further simplify our 
functions by considering the two regimes of interest we observe in our models: one where 
we will have </> = (/>o today if we intend to recover G e s = G at the present; and another one 
where we do not recover <f>g at the present, meaning that, essentially, we instead recover the 
matter domination attractor solution for <p at early times given by Eq. (2.9). For the first 
case, we obtain: 


fQSi («) 


8 - 6Q m 3\/(3y 7 d(l — fl m ) (1 + w e s) 3 ^ 2 

d{ 2-n m ) 3d (1 + w e g) — 2 


while for the second case we find 


£qS 2 ( a ) 



V6y/d(l + w eB ) 3 / 2 
(-2+3d(l+™ eff ))m eff 

(«) 


(5.10) 


(5.11) 


The GR limit of £qs = 1 is recovered in both situations when we take wbd —> oo. For the 
first case, as expected, £qs = 1 when a = 1 since we have <f> = 4>g today. In the second 
case, £qs < 1 today since the present-day value of the scalar field in these circumstances 
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Figure 7 : We show the late-time evolution of G e g /G in the quasi-static regime using Eq. (4.9) 
in the upper plot. The solid lines show G e g/G with 4> = 4 >o today, while the dashed lines 
do not have that restriction. The red lines are for <ubd = 100, the blue lines for cjbd = 500 
and the green lines have wbd = 1000. The bottom plot shows the relative error of our 
approximations to the exact numerical solutions in %. 


will always be larger than fo. This can be observed in Fig. 7, where we compare these 
approximations to the exact numerical solution of £ and we see they work considerably well. 

We can now try and understand the dependence of £qs on the different parameters. 
Looking at d = 2u>bd + 3, it becomes clear that increasing the Brans-Dicke parameter leads 
1 ° £qs becoming closer to 1 throughout the late-time cosmological evolution: its slope at 
a = 1, as given by Eq. (5.10), decreases since it depends on the inverse of d or \fd. Then, 
looking at Eq. (5.11), we see that the present-day value of £qs increases towards 1 due to 
the exponents of the terms shown becoming extremely small. 

Looking at the dependence of £qs on fi m , we realize it is similar to that on wbd- Increas¬ 
ing Q m leads to both the slope of £qs decreasing in Eq. (5.10) as well as the present day-value 
tending to 1 in Eq. (5.11). In Eq. (5.11) we also see that, remarkably, our approximation 
recovers the matter dominated attractor solution value of £qs = 1 when —> 1. 

Lastly, we have the effective equation of state parameter, w e fj. Looking at Eq. (5.10), 
we see that the slope of £qs will increase as w e s becomes less negative, making its evolution 
more noticeable for larger w e g- when all other parameters remain fixed. Also, the exponent of 
the second term in Eq. (5.11) increases for large w e g, leading to a significant departure of £qs 
today, producing values of £qs(ci = 1 ) that are detectably smaller than 1 in a clear departure 
from standard GR. This is a reflection of the effect of increasing rc e ff on the evolution of the 
scalar field </>: the higher w e g- is, the sooner cj) departs from the matter domination attractor 
solution and the larger its present-day value will be. 

In Fig. 8, we plot £qs as a function of wbd at a = 1, using Eq. (5.11). We see that if we 
don’t fix cj) = (f>o today, there is a significant, possibly detectable, deviation from the standard 
GR value, £qs = 1, even for very large wbd- Of course, we also see that when wbd —> oo, £qs 
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Figure 8: We show the evolution of £qs at a = 1 for S7 m = 0.308 and w e g = — 1 as a function 
of wbd- For this plot we have used Eq. (5.11), therefore assuming that </>(a = 1) may not be 
equal to in¬ 


tends to 1. Therefore, in order to be competitive with Solar System constraints cubd > 10 4 
[18, 19], we would have to able to measure £qs with a precision of around 10 -4 . 

6 Discussion 

In this paper we have applied the designer approach to the extended Brans-Dicke theory with 
the explicit presence of a self-interacting potential V ( q i>). By fixing the expansion history to 
that of an effective rcCDM dark energy model, we are able to retrieve the scalar field evolution 
under the assumption that the main contribution to the effective dark energy density comes 
from the potential V(4>). 

The numerical solutions we obtain have the property of respecting the matter domina¬ 
tion attractor solution of Brans-Dicke models at early-times. At late-times, the scalar field 
departs from this solution and evolves more rapidly and towards larger values, yielding a 
value today larger than <po, where is the present-day value of the matter regime attractor 
solution that ensures that one would measure G e fj today equal to the actual Newton’s grav¬ 
itational constant, G, in a matter-dominated Universe. This transition from the attractor 
solution happens earlier whenever we take a larger dark energy equation of state, u> e ff. How¬ 
ever, if we constrain the present-day value of 4> to be equal to (f> o, our numerical solutions 
follow the power-law behavior of the attractor solution, shifted towards smaller values at 
early-times. When the evolution departs from the matter dominated behavior, we are then 
able to recover </>(a = 1) = </>o as intended. 

We were able to obtain separate analytical approximations for the evolution of the 
scalar field when w e g = — 1 and w e g > — 1, which we then used to construct a global solution 
valid for w e g > — 1. These approximations work remarkably well, with errors of sub-percent 
for large values of u>bd- These approximations also allowed us to reconstruct the late-time 
functional form of the potential U(</>); we found a simple run-away potential whose slope 
is inevitably dependent on u; e ff and cjbd- We reiterate that we have limited our analysis 
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to constant w e f[ so as to obtain analytical solutions which will shed light on the parameter 
dependence of the various observables we are considering; a non-constant w e g will severely 
complicate any attempts at doing so. However, we stress that the numerical implementation 
of the designer approach presented in Sec. 3 can be easily extended to a non-constant u> e ff. 

With these analytic approximations in hand, we then focused on the phenomenological 
parameters that describe the sub-horizon evolution of the linear perturbations of the theory. 
We showed how the effective scale of the theory, which we designated by &m, is of order 
the cosmological horizon; as a result we find that there is negligible scale dependence of the 
phenomenological parameters on observable scales. We found that the ratio between the 
Newtonian potentials, r/ = T /<3? is constant throughout the cosmological evolution, for large 
values of the Brans-Dicke parameter [37]. We also found simple analytical expression for 
£ = G e ff/G which depend explicitly on the parameters of the theory, as seen in Eqs. (5.10) 
and (5.11). 

One of the main features of this model is the possibility of having £qs / 1 today; this 
is due to the departure of the scalar field from the matter dominated attractor solution at 
late-times such that its present-day value will be larger than </>o- The present-day value of 
£qs at a = 1, given in Eq. (5.11) tends to 1 as u>bd —> oo since the exponent of the terms 
shown tend to zero. Also, as for w e ff > —1, the exponent of one of the terms increases, 
leading to smaller values of £qs today, even when cjbd is very large. 

If, however, we impose £qs to be 1 today, its evolution is predicted by Eq. (5.10). In 
these circumstances, the main distinguishing point between this model and standard GR will 
be the slope of £qs at the present: for the extended Brans-Dicke theory it can be different 
from zero. We note that, even when rc e ff = — 1, the predicted slope is different from zero. 
Hence, even the simple extended Brans-Dicke+ACDM model could be ruled out if £qs is 
found to not vary close to the present. 

Finally, we note that in order to attain constraints on cjbd that are competitive with 
those obtained in Solar-system tests [18, 19], £qs and t]qs would naively need to be con¬ 
strained with a precision of around ICE 4 . This is a formidable challenge, but one should 
bear in mind that t]qs ^ 1 throughout (at least) the matter dominated era while the same is 
possible for £,qs- This means that there will be a cumulative effective (as shown in [39, 40]) 
which means that constraints on the growth rate (or weak lensing) of order 10 -3 or even 
ICE 2 might be sufficient to place competitive constraints on tcbd- 
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A Correction factor for V(0) 


In this appendix we show the correction factor one can add to the potential V(<j>) in order to balance 
the scalar field dynamics in the exact numerical evolution of in order to recover a flat Universe 
today. Effectively, we want to solve the equation 


<t>'(ao) cubd /£M\ 2 1 - .Dff m _ 

</>(a 0 ) 6 V </H a o) J <Ma 0 ) 


(A.l) 


where ao = 1 and D will be the correction factor such that 1 — Dtt m = 12^, as discussed in Sec. 3. 
First we show that factor using our solution for w e s = — 1 using Eq. (4.2): 


D = 


H a i)g(a 0 ) 


3 d (fi m - 2) (8 - 6f2 m + d (2 - n m ) (1 - Sl m )) + 2wbd (4 - 3fi m )‘ 


n„ 


g(ai) 


3 n m d 2 (fi m - 2 )- 


(A.2) 


where g(a) is defined in Eq. (4.3) and 4>{ai) is the value of the scalar field at the starting redshift o,. 

Lastly, we show the correction factor for the case w e s > — 1. For this part we have used the 
late-time solution for </> given by Eq. (4.6): 


1 4>{a,i) f 0 3vW dsj 1 - U m (1 + w eS ) 3/2 9w B d (1 - f^m) d(l + w e ff) 3 \ M 0 s 

_ On L! m m+ — 2 + 3d(1 + w e s) (2 — 3c?(1 + w e ff)) 2 ) k 

where c = /(a 0 )//(ai), and /(a) is defined by Eq. (4.7). 
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